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Abstract 

We have developed a scaled parametric equation of state to describe and 
predict thermodynamic properties of supercooled water. The equation of 
state, built on the growing evidence that the critical point of supercooled 
liquid-liquid water separation exists, is universal in terms of theoretical scaling 
fields and is shown to belong to the Ising-model class of universality. The 
theoretical scaling fields are postulated to be analytical combinations of the 
physical fields, pressure and temperature. The equation of state enables us 
to accurately locate the "Widom line" (the locus of stability minima) and 
determine that the critical pressure is considerably lower than predicted by 
computer simulations. 

Upon supercooling, water exhibits anomalous behavior with sharply increasing heat ca- 
pacity, isothermal compressibility, and the magnitude of negative thermal expansivity [1]. A 
thermodynamically consistent view on the global phase behavior of supercooled water was 
formulated by Poole et al. [2] . According to this view, the observed anomalies are associated 
with density and entropy fluctuations diverging at a critical point of liquid-liquid coexistence 
that terminates the line of flrst-order transitions between two liquid aqueous phases: low- 
density liquid and high-density liquid. This " second-critical-point" scenario is supported by 



extensive Monte Carlo and molecular dynamic simulations [3] , by a modified van der Waals 
model that includes hydrogen-bond contributions [4] , and by the limited but impressive ex- 
perimental evidence [5,6]. An alternative interpretation of the phase behavior of water, the 
"singularity-free" scenario, attributes the increase in response functions upon supercooling, 
through a thermodynamic consistency argument, to the existence of a negatively-sloped 
locus of density maxima in the P-T plane [7] . However, the second-critical-point scenario 
seems more plausible in view of the experimental evidence of a first-order phase transition 
between two amorphous- water glasses [1,5,6]. The global phase diagram of water is shown 
in Fig. 1. Two remarkable features make the second critical point in water phenomeno- 
logically different from the well-known gas-liquid critical point. The negative slope of the 
liquid-liquid phase transition line in the P-T plane means that the higher density liquid 
water is the phase with larger entropy. A very large value of the slope at the critical point 
(about 30 times larger than for the vapor-liquid transition at the critical point) indicates the 
significance of the entropy change with respect to the density change, and correspondingly, 
the importance of the entropy fluctuations. However, the location of the liquid-liquid critical 
point, especially the value of the critical pressure is uncertain. The simulation data yields a 
variety of the critical-pressure values, from negative pressures to 3.4 kbar [1,8]. 

In this Letter, we present a scaled parametric equation of state to describe and predict 
thermodynamic properties of supercooled water. The equation of state is universal and 
belongs to the Ising-model class of universality. The equation of state enables us to accu- 
rately locate the critical point and the "Widom line" [9] (the locus of stability minima and 
order-parameter fluctuation maxima). In particular, we conclude that the critical pressure 
is considerably lower than obtained by computer simulations and we predict thermodynamic 
properties in the regions inaccessible to experiments. 

It is commonly accepted that the critical behavior of all fluids, simple and complex, 
belong to the universality class of the three-dimensional Ising model [10]. Water is not 
an exception [11]. Near the critical point the critical (fluctuation-induced) part, "^^r, of 
an appropriate fleld-dependent thermodynamic potential ^ is a universal function of two 
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scaling fields, "ordering", hi, and "thermal", /i2 [10]: 



where a — 0.109, /3 = 0.326, 7 = 1.239, are universal exponents (interrelated as a + 
2/9 + 7 = 2) in the scaling power laws (as a function of /i2 at hi = 0) for the "weak" 
susceptibility, the order parameter and the "strong" susceptibility, respectively. The first 
derivatives of the thermodynamic potential with respect to the scaling fields define two 
scaling densities, the "order parameter" 0i = —d'^cr/dhi and the "thermal density" 02 = 
—d'^cr/dh2. Respectively, the second derivatives define three susceptibilities, "strong" 
Xi = id(f)i/dhi)f^^, "weak" X2 = (902/9/12)^^, and "cross" X12 = id(f)i/dh2)h^ = (902/9/ii)^^. 
The universal scaling function / contains two system-dependent amplitudes that originate 
from the initial Hamiltonian. In the mean- field approximation [a — 0, /3 — 1/2, 7=1), 
the critical part of the thermodynamic potential is represented by a Landau expansion, 
l'^o0i ~ hicpi, where the constants ag and Uq (the same coefficients as in 
the initial Hamiltonian) play the role of the mean-field system-dependent amplitudes. 

To apply the universal expression (1) for describing the liquid-liquid critical point in 
supercooled water, we assume the scaling fields are analytical combinations of physical 
fields, the pressure P and the temperature T: 

hi = aiAP + aaAT + agAP^ (2) 
h2 = 61 AT + 62AP (3) 

with AP = (P - Pc) / (pcRTc) and AT = (T - %) (where the subscript "c" here and 
below indicates the critical parameters) and and hi are system-dependent coefficients. 
Representation of scaling fields through linear mixing of physical fields is commonly used to 
incorporate asymmetric fiuid criticality into the symmetric Ising model [?,13]. To account 
for the strong curvature of the liquid-liquid transition line, defined as hi = 0, we added a 
non- linear pressure term, 03 AP^, in hi. 
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Such a representation of the scahng fields corresponds to the practically convenient choice 
of the molar Gibbs energy (chemical potential) as the field-dependent potential so that 
^ = ^(P, T) and ^'cr = ^cr(/ii, /i2)- In this formulation we neglect the difference be- 
tween —Ay = — {V — Vc) /Vc — {p — Pc) / p and {p — pc) / Pc, a reasonable approximation 
for weakly compressible liquids. Any two coefficients in the scaling fields may be absorbed 
in two system-dependent amplitudes of the scaling function /, so that we adopt 02 = 1 and 
62 = —1- The negative sign of 62 indicates that the liquid-liquid phase separation in super- 
cooled water occurs with increase of pressure (Fig. 1), in contrast to the vapor- liquid phase 
separation. The value of ai, a^, and bi can be determined from the shape of the liquid- 
liquid first-order transition curve. The "Widom line" in the one-phase region, hi — 0, is an 
analytical continuation of the liquid- liquid transition curve from C to lower pressures and 
higher temperatures. 

The mixing of physical fields into the scaling fields, defined by Eqs. (2) and (3), means 
that the order parameter is a combination of molar entropy and molar volume. In the linear 
approximation 

61AV' + h2AS 



0262 + (ai)efr ^1 ' 



■>2 = 



(4) 
(5) 



0262 - (ai)eff^l' 

where A^ = {S — Sc)/R, with R being the gas constant, and (ai)eif = [dh\/dP^ _ = 
oi 2a3AP. 

As far as the physical fields are mixed into the scaling fields, the physical properties, 
such as the isobaric heat capacity Cp, the isothermal compressibility kt, and the thermal 
expansivity ctp, will not exhibit universal power laws when measured along isotherms or 
isobars; instead, their apparent behavior will be determined by a thermodynamic path and 
by the values of the mixing coefficients in Eqs. (2) and (3). As follows from Eqs. (2) and (3), 
the critical (fiuctuation induced) parts of the dimensionless isobaric heat capacity, isothermal 
compressibility, and thermal expansivity are expressed through the scaling susceptibilities: 
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= -T {alxi + 2a2&iXi2 + &?X2) , 
1 / dV 
V \ dP 



(^^)cr -ikT\ (7) 



T 



= i ((«l)efr Xl + 2 (ai)eff ^'2X12 + &2X2) , 
(«^)cr=ife) -(«p)b (8) 



V \ dT 



p 



1 

= ((«i)eff «2Xi + ((ai)eff ^1 + 02^2) X12 + bib2X2) ■ 

where T — T/Tc, P — P/pcRTc, and the subscript "b" indicates the property backgrounds. 

We use the simplest form of a scaled parametric equation of state, the so-called "linear 
model", which represents the scaling fields and scaling susceptibilities as functions of the 
"polar" variables r and 6 [13,14]: 

hi = ar^^+W (1 -e^), h2 = r{l- hH^) , (9) 

Xl = -r-^ci(^), X12 = A;r^"Vi2(^), X2 = akr-''c2{e) - B,, 
a 

where the coefficient 6^ = (7 — 2/3) /^(l — P) ~ 1.36 is a universal constant, while a and 
k are system-dependent amplitudes, and B^r is the so-called "critical background" of order 
ak [13]. The analytical functions Ci{6), 02(0), and CuiO) are calculated in ref. [13]. A 
remarkable feature of the "linear model" is that the singularities in the thermodynamic 
functions are only related to the variable r, while the properties are analytical with respect 

to e. 

This model offers a consistent scaling description of the available experimental data in 
supercooled water. Using high-resolution experimental heat-capacity data [15] shown in 
Fig. 2, we optimized the location of the critical point and the system dependent ampli- 
tudes. Based on the most recent estimate of the liquid-liquid phase transition curve given 
by Mishima [6], we have obtained the coefficients ai = bi = 0.0078 and 03 = 0.062. These 
particular numbers correspond to Pc — 27 MPa, the value optimized by our equation of state. 



We assume that the "Widom hne" is described by the same coefficients. Furthermore, to 
reduce the number of adjustable parameters, we assume that the ratio fc/a = 1, as obtained 
for the three dimensional Ising model with short range interactions [16]. Hence, only two 
adjustable parameters, namely, Pc and a — k, have been used to describe the anomalous 
parts of the thermodynamic properties. The non-critical background of the heat capacity 
was approximated as a linear function of temperature. We obtained a = k = 0.47 and 
Pc — 27 MPa with the critical temperature corresponding to this pressure Tc — 232 K. The 
critical point, obtained from our equation of state, is located at a much lower pressure than 
previously predicted from computer simulations (see Fig. la). 

With the given amplitudes and location of the critical point, we predict the behavior of 
the compressibility and expansivity, shown in Figs. 3 and 4, by adjusting only their non- 
critical backgrounds. The molar volume as a function of temperature was taken from ref. 
[17]. The predictions appear to have excellent agreement with the experimental data [18,19]. 
While it is difficult to establish the error bars for the obtained Pc value, the parametric 
equation of state certainly excludes the critical pressure above 50 MPa or below 10 MPa. 
We also conclude that the mean-field scenario is unlikely. The mean-field scenario cannot 
predict the anomalous behavior of isothermal compressibility within our model. While the 
major contribution in the heat-capacity anomaly is strong susceptibility, xi, (6i is small) 
diverging both in mean-field and in scaling theory, the major contribution in the isothermal 
compressibility anomaly is the weak susceptibility X2, («i is small) which shows no anomaly 
in mean-field approximation. The major contribution in the critical part of the expansivity 
comes from the cross susceptibility Xi2 as both ai and bi are small. These features make 
the second critical point in water essentially different from the liquid-vapor critical point 
where Cp, kt, and ap all diverge strongly, as Xi, and from the liquid-liquid critical points 
in binary fluids where Cp, kt, and ap all diverge weakly, as X2- 

There are obvious limitations of our equation of state. First, the model used in this 
work is accurate only asymptotically close to the critical point (r << 1) while all mea- 
surements in supercooled water have been taken far beyond the asymptotic region. The 
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experimental range of r, the parametric distance to the critical point, may be as large as 0.5. 
However, this is the first estimate of the critical parameters for the second critical point in 
water based on experimental data, and not on computer simulations of "water like" models. 
Including non-asymptotic corrections to the parametric equation of state would change the 
adjustable backgrounds while not significantly affecting the critical parameters. To more 
accurately describe and predict the properties in a broader range of pressures and densities 
in supercooled water, a "global" crossover equation of state [14], based on a rehable mean- 
field equation of state, such as a modified van der Waals model [4], is required. Moreover, 
we did not address an intriguing possibility of the existence of multiple critical points in 
supercooled water, as predicted by some simulated water models [20]. 

In this work, the order parameter is phenomenologically expressed through molar volume 
and entropy, with entropy being the major contribution. A clarification of the relation 
between this phenomenology and the microscopic nature of the order parameter [21] would 
help in better understanding the physics of phase transitions in supercooled water. 

We acknowledge valuable discussions with P. Debenedetti, O. Mishima, H.E. Stanley, 
and B. Widom. 



7 



REFERENCES 



[1] P.F. Debenedetti, J. Phys.: Condens. Matter 15, R1669 (2003). 

[2] P.H. Poole, F. Sciortino, U.Essmann, and H.E. Stanley, Nature 360, 324 (1992). 

[3] S. Sastry, F. Sciortino, and H.E. Stanley, J. Chem. Phys. 98, 9863 (1993); P.H. Poole, F. 
Sciortino, U. Essmann, and H.E. Stanley, Pys. Rev. E 48, 3799 (1993); P.H. Poole, U. 
Essmann, F. Sciortino, and H.E. Stanley, Phys. Rev. E 48, 4605 (1993); H.E. Stanley, 
C.A. Angell, U. Essmann, M. Hemmati, P.H. Poole, and F. Sciortino, Physica A 205, 
122 (1994). 

[4] P.H. Poole, F. Sciortino, T. Grande, H.E. Stanley, and C.A. Angell, Phys. Rev. Lett. 
73 1632 (1994). 

[5] O. Mishima and H.E. Stanely, Nature 392, 164 (1998). 

[6] O. Mishima, J. Chem. Phys. 23 154506 (2005); O. Mishima (personal communication). 

[7] S. Sastry, P.G. Debenedetti, F. Sciortino, and H.E. Stanley, Phys. Rev. E 53, 6144 
(1996); P.G. Debenedetti, Nature 392, 127 (1998); L.P.N. Rebelo, P.G. Debenedetti, 
and S. Sastry, J. Phys. Chem. 109, 626 (1998). 

[8] S. Harrington, R. Zhang, P.H. Poole, F. Sciortino, and H.E. Stanley, Phys. Rev. Lett. 
78, 2409 (1997); F. Sciortino, P.H. Poole, U. Essmann, and H.E. Stanley, Phys. Rev. E 
55, 727 (1997); M. Yamada, S. Mossa, H.E. Stanley, and F. Sciortino, Phys. Rev. Lett. 
88, 195701 (2002); A. Scala, F.W. Starr, E. La Nave, H.E. Stanley, and F. Sciortino, 
Phys. Rev. E 62, 8016 (2000); H. Tanaka, Nature 380, 328 (1996); H. Tanaka, J. Chem. 
Phys. 105, 5099 (1996). 

[9] L. Xu, P. Kumar, S.V. Buldyrev, S.-H. Chen, P.H. Poole, F. Sciortino, and H.E. Stanley, 
Proc. Natl. Acad. Sci. 102, 16558 (2005). 

[10] M.E. Fisher in Critical Phenomena, edited by F.J.W. Hahne, Vol. 186 (Springer, Berhn, 



8 



1982), p. 1. 

[11] M.A. Anisimov, J.V. Sengers, and J.M.H. Levelt Sengers, in Near-Critical Behavior of 
Aqueous Systems, edited by D.A. Palmer, R. Fernandez-Prini, A.H. Harvey, The Phys- 
ical Properties of Aqueous Systems at Elevated Temperatures and Pressures: Water, 
Steam and Hydrothermal Solutions (Academic Press, 2004). 

[12] N.B. Wilding and A.D. Bruce, J. Phys.: Condensed Matter 4, 3087 (1992); M.E. Fisher 
and G. Orkoulas, Phys. Rev. Lett. 85, 696 (2000); M.A. Anisimov and J.T. Wang, Phys. 
Rev. Lett, (in press). 

[13] M.A. Anisimov, V.A. Agayan, and P.J. CoUings, Phys. Rev. E, 57, 582 (1998). 

[14] M.A. Anisimov and J.V. Sengers, in Equations of State for Fluids and Fluid Mixtures, 
edited by J.V. Sengers, R.F. Kayser, C.J. Peters, and J.J. White Jr. (Elsevier, Amster- 
dam, 2000), p. 381. 

[15] C.A. Angel, W.J. Sichina, and M. Oguni, J. Phys. Chem. 86, 998 (1982). 

[16] Y.C. Kim, M.A. Anisimov, J.V. Sengers, and E. Luijten, J. Stat. Phys. 110, 591 (2003). 

[17] B.V. Zheleznyi, Russ. J. Phys. Chem. 43, 1311 (1969). 

[18] H. Kanno and C.A. Angel, J. Chem. Phys. 70, 4008 (1979). 

[19] D.E. Hare and CM. Sorensen, J. Chem. Phys. 84, 5085 (1986). 

[20] V.B. Henriques, N. Guisoni, M.A. Barbosa, M. Thielo, M.C. Barbosa, Mol. Phys. 103, 
3001 (2005); I. Brovchenko, A. Geiger, A. Oleinikova, J. Chem. Phys. 123, 044515/1 
(2005). 

[21] T.M. Truskett, P.G. Debenedetti, S. Sastry, S. Torquato, J. Chem. Phys. Ill 2647 
(1999). 



9 



Captions of figures 

Figure 1. (a) Phase diagram for water with vapor-hquid and hquid-hquid critical points. 
Solid curves are vapor-liquid, hquid-liquid, and liquid-solid phase transitions; dashed is the 
Widom line; dotted is the compressibility maxima; crosses are literature estimates for the 
liquid-liquid critical point. C is the vapor-liquid critical point; C is the liquid-liquid critical 
point predicted by our equation of state, (b) The vicinity of the liquid- liquid critical point. 
[5,8]. 

Figure 2. Heat capacity measurements (stars) [15] compared to the heat capacity pre- 
dicted by our scaling parametric equation of state (solid curve) for a critical point of Pc — 27 
MPa, Tc = 232 K. Non-critical background plotted as a thin solid curve. 

Figure 3. Isothermal compressibility experimental data [18] at 10 MPa (o), 50 MPa 
(A), 100 MPa (O), 150 MPa (*), and 190 MPa (□), compared with our prediction at the 
same pressures (solid curves). 

Figure 4. Thermal expansivity experimental data [19] (a), at ambient pressure, com- 
pared with our prediction (solid curve). Non-critical background plotted as a thin solid 
curve. 
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